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Abstract 


Experimental  evidence  regarding  pressure  gradients  in  Bridgman  anvil  devices 
is  reviewed  and  some  new  experimental  data  are  presented.  Although  there  are 
many  differences  in  details,  all  the  experimental  evidence  is  found  to  fit  a  general 
pattern.  A  simple  empirical  relationship  describing  the  pressure  distribution  is 
derived  which  is  consistent  with  experimental  observation.  The  experimentally 
determined  pressure  is  used  to  estimate  errors  introduced  into  various  types  of 
anvil  experiments  by  pressure  gradients.  Numerical  values  for  the  correction 
factors  are  given  in  graphical  form  for  a  wide  range  of  experimental  conditions. 
Analysis  of  shear  experiments  indicates  that  pressure  gradients  are  much  smaller 
under  dynamic  conditions.  Analysis  also  shows  that  the  peripheral  sample  region 
which  is  partially  non- load- supporting  frequently  has  a  greater  detrimental  effect 
on  anvil  experiments  than  does  the  pressure  gradient  which  exists  over  the  central 
portion  of  the  sample  cell. 
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Pressure  Gradients  in  Bridgman  Anvil  Devices 


I.  INTRODUCTION 


Opposed  anvils  are  the  simplest  and  most  widely  used  form  of  very-high-pres- 
sure  apparatus.  Bridgman  (1935,  1936,  1937)  used  them  extensively  in  the  vicinity 
of  room  temperature  and  they  are  now  being  used  at  high  and  low  temperatures  as 
well.  They  have  been  used  to  investigate  a  wide  range  of  physical  phenomena. 
Although  it  has  long  been  realized  that  truly  hydrostatic  pressures  were  lot  ob¬ 
tained  with  anvils,  the  magnitude  of  the  pressure  gradients  which  can  exist  in  this 
type  of  apparatus  was  not  appreciated  fully  until  recently.  The  nature  of  the  pos¬ 
sible  pressure  gradients  and  their  influence  on  certain  types  of  measurements 
constitute  the  subject  of  this  paper. 

For  purposes  of  discussion  it  is  convenient  to  divide  opposed  anvil  devices 
into  categories  -  static  and  dynamic.  Most  opposed  anvil  devices  are  static,  which 
means  that  the  anvils  are  not  rotated  relative  to  each  other.  A  relatively  small 
number  are  dynamic,  in  which  one  anvil  is  rotated  or  oscillated  while  the 
other  is  held  stationary.  The  dynamic  type  is  used  primarily  for  shear  strength 
measurements  or  in  phase  equilibrium  investigations. 

The  pressure  distribution  in  a  sample  compressed  between  anvils  may  depend 
significantly  on  whether  the  anvils  are  static  or  dynamic.  It  will  also  depend  on 
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the  nature  of  the  sample,  the  anvil  material,  the  temperature  and  pressure  applied, 
the  geometrical  arrangement,  and  other  f  ._tors.  In  a  dynamic  device  the  pressure 
distribution  may  depend  on  the  rate  of  strain  of  the  sample  material. 

The  composition  and  structure  of  the  wafer-like  sample  cell  used  in  an  anvil 
apparatus  vary  markedly  from  one  experiment  to  the  next.  In  shear  strength 
measurements  and  some  geocnemical  investigations  the  sample  cell  consists  of  a 
homogeneous  disc,  but  this  simple  configuration  is  the  exception  rather  than  the 
rule.  Most  experimental  arrangements  consist  of  a  central  disc  surrounded  by 
an  annular  ring  which  retains  the  sample  between  the  anvils.  Pyrophyliite  is 
commonly  used  for  the  retaining  ring,  although  stainless  steel  or  nickel  are  some¬ 
times  employed.  The  properties  and  dimensions  of  the  retaining  ring  relative  to 
those  of  the  central  disc  have  an  important  bearing  on  the  pressure  distribution. 

The  central  disc  in  the  sample  cell  may  be  the  material  under  investigation 
or  it  may  serve  only  as  a  pressure  transmitting  medium.  In  the  latter  case  the 
sample  itself  may  be  a  wire  or  foil  on  which  electrical  or  other  measurements  are 
to  be  made.  The  presence  of  such  a  specimen  embedded  in  the  central  disc  further 
complicates  the  pressure  distribution.  This  marked  variation  in  the  composition 
of  the  sample  cell  greatly  complicates  the  problem  of  determining  the  pressure 
gradient  between  opposed  anvils,  and  no  doubt  is  partly  responsible  for  the  widely 
different  gradients  reported. 

Bridgman  (1937)  calculated  the  deformation  of  anvils  under  pressure  using 
classical  elasticity  theory  assuming  a  uniform  pressure  distribution. 

His  analysis  predicted  that  the  anvils  and  sample  wafer  would  assume  a  len¬ 
ticular  shape  when  compressed.  Because  the  permanent  deformation  of  anvils 
and  sample  wafers  conformed  to  his  prediction,  he  assumed  that  a  uniform  pres¬ 
sure  distribution  was  a  atisfactory  first  approximation.  Bridgman  also  believed 
that  a  more  careful  analysis  would  show  that  the  maximum  pressures  occurred  in 
an  annular  region  concentric  with  the  axis  of  compression.  Recent  experimental 
work  has  shewn  this  qualitative  description  to  be  valid  in  some  instances  but  not 
in  others. 

Christiansen  et  al  (1962)  determined  the  pressure  distribution  between  anvils 
using  silica  glass  as  a  pressure  indicator.  They  established  that,  on  compression 
at  room  temperature,  this  material  experiences  a  permanent  densification  which 
is  linearly  related  to  the  pressure.  In  compression  tests  on  0.  24-in.  -diam  glass 
wafers,  the  pressure  near  the  edge  of  the  disc  was  much  larger  than  at  the  center, 
100  kilobars  vs  70  kilobars  in  a  typical  case  when  the  nominal  applied  pressure  was 
about  73  kilobars.  In  another  series  of  tests,  in  which  the  glass  disc  was  encircled 
by  a  pipestone  retaining  ring,  the  average  pressure  on  the  glass  disc  was  greater 
than  the  nominal  applied  pressure  by  from  10  to  25  percent.  Unfortunately  the 
pressure  gradient  was  not  determined  in  the  latter  experiments. 
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Although  these  results  may  be  considered  typical  because  they  demonstrate 
the  occurrence  of  substantial  pressure  gradients  and  also  tend  to  confirm  Bridgman’s 
qualitative  description,  they  are  unusual.  In  most  experimental  configurations, 
the  maximum  pressure  occurs  at  the  center  of  the  anvils.  Pressure  multiplication 
ei/ects  have  been  reported  in  which  the  pressure  at  the  center  is  as  much  as  2  or 
2.  5  times  the  nominal  ipplied  pressure. 

For  example,  Myers  et  al  (1963)  made  a  rather  extensive  investigation  on  a 
sample  wafer  design  consisting  of  a  pyrophyliite  disc  surrounded  by  a  nickel  re¬ 
taining  ring.  Small  calibrant  wires  were  located  at  various  radial  positions.  Tne 
wi^es  were  fabricated  from  materials  having  known  phase  transitions  which  can 
be  detected  electrically.  Thus  the  nominal  applied  pressure  required  to  produce 
a  known  transition  pressure  could  be  measured  as  a  function  ot  xne  radial  position 
of  the  caHbrant  wire.  The  measurements  did  not  determine  the  pressure  gradient 
directly;  however  they  (lid  demonstrate  that  substantial  pressure  gradients  must 
exist.  In  particular  it  was  found  that  the  pressure  at  the  center  of  the  anvils  might 
be  either  greater  or  less  than  the  nominal  applied  pressure,  depending  on  the 
diameter-to-thickness  ratio  of  the  sample  wafer.  For  thick  samples  a  pressure 
multiplication  of  as  much  as  2.  5  was  observed  at  the  center  of  the  anvils  while  for 
thin  samples  a  pressure  deficiency  was  noted.  The  multiplication  effect  was  ob¬ 
served  over  a  wide  range  of  pressures,  and  the  ratio  of  nominal  applied  pressure 
to  actual  pressure  at  the  center  of  the  anvils  was  practically  constant  in  the  range 
up  to  90  kilobars. 

Montgomery  et  al  (1963)  have  investigated  the  pressure  gradient  in  a  silver 
chloride  disc  surrounded  by  a  pyrophyliite  retaining  ring.  They  report  a  depen¬ 
dence  of  the  central  pressure  on  the  diameter-to-thickness  ratio  which  is  quali¬ 
tatively  the  same  as  that  found  by  Myers  et  al.  Using  calibrant  wires  located  at 
various  radial  positions  they  also  demonstrated  that  a  large  pressure  gradient 
existed  in  the  silver  chloride  discs.  The  nominal  applied  pressure  required  to 
produce  the  bismuth  I- II  transition  was  related  lineaily  to  the  radial  position  of 
the  calibrant  wire.  Similar  observations  developed  using  the  bismuth  VI-VUI  tran¬ 
sition.  The  fit  to  the  straight  line  is  very  good,  the  experimental  reproducibility 
being  a  few  percent.  It  is  important  to  note  that  this  does  iK>t  necessarily  imply 
a  linear  pressure  gradient  for  reasons  which  will  be  discussed. 

Montgomery  et  al  also  showed  that  manganin  wires  in  the  form  of  circular  arcs 
concentric  with  the  center  of  the  anvils  gave  very  reproducfole  plots  of  relative 
resistance  vs  nominal  applied  pressure.  Hysteresis  observed  remained  within  the 
uncertainty  of  the  measurements,  and  wires  of  different  radii  and  arc  length  yielded 
resistance  plots  of  similar  shape.  Since  manganin  is  very  strain- sensitive,  these 
results  showed  clearly  that  the  pressure  distribution  had  the  avial  symmetry  ex¬ 
pected. 
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It  was  also  discovered,  by  making  a  single  plausible  assumption  that  the  pres¬ 
sure  at  any  radial  position  was  always  proportional  to  the  nominal  applied  pres¬ 
sure,  that  all  relative  resistance  data  obtained  with  wires  of  different  radii  could 
be  superimposed  to  form  a  single  curve  of  relative  resistance  versus  pressure. 

The  assumption  was  confirmed  for  the  center  of  the  anvils  by  Myers  et  al,  as  noted 
previously.  The  calibration  curve  was  used  to  calculate  the  pressure  distribution 
between  the  anvils  at  several  pressures  and  it  was  found  that  the  central  pressure 
was  less  than  the  nominal  applied  pressure  and  that  the  pressure  increased  lin¬ 
early  with  radial  position. 

In  orde^  to  make  the  calculations,  an  additional  plausible  assumption  was  made, 
that  the  pressure  dropped  abruptly  in  a  linear  manner  to  some  snii.ll  value  com¬ 
parable  to  the  shear  strength  cl  the  retaining  ring  in  the  region  near  the  periphery 
of  the  sample.  The  width  of  this  region  was  about  0.  030  in.  Thns  on  1/2-in.  - 
diam  anvil e  an  annular  ring  comprising  about  23  percent  of  the  total  anvil  area 
supported  only  a  small  fraction  of  the  applied  load.  Also,  the  size  of  this 
transition  region  was  independent  of  the  anvil  size  and  the  applied  load. 

The  existence  of  this  border  region  is  also  evident  from  a  series  of  compres¬ 
sion  tests  on  homogeneous  wafers  of  several  materials  made  by  Vaisnys  and 
Montgomery  (1964).  They  discovered  that  the  load-bearing  area  of  the  anvils  varied 
from  100  percent  to  50  percent  of  the  geometrical  area,  depending  on  the  material 
compressed  and  on  the  diameter-to-thickness  ratio.  Thick  samples  resulted  in 
the  smallest  load-bearing  areas  in  agreement  with  the  pressure- multiplication 
effect  observed  by  others. 

The  pressure  distribution  between  miniature  diamond  anvils  has  been  observed 
optically  by  Duecker  and  Lippincott  (ASME  Paper,  to  be  published)  on  homogeneous 
discs.  A  pressure  multiplication  effect  of  two  or  more  was  found  at  low  pressures 
but  the  effect  diminished  at  high  pressures.  However,  the  constant  ratio  of  central 
pressure  to  nominal  applied  pressure  observed  by  Myers  et  al  in  the  high-pressure 
region,  and  at  sumed  by  Montgomery  et  al,  was  not  observed.  At  low  loads  the 
radial  pressure  distribution  appeared  approximately  parabolic  and  a  theoretical 
model  was  given  which  correlated  well  with  the  observations.  At  nominal  applied 
pressures  of  20  kilobars  or  more,  the  observations  fit  a  linear  pressure  distribution 
rather  well  with  the  exception  oi  a  small  area  at  the  certer  of  the  anvils.  Since 
the  pressures  near  the  anvil  edges  were  on  the  order  of  a  kilobar,  the  observed 
pressure  gradients  were  gigantic. 

It  remains  difficult  to  draw  any  general  conclusions  from  the  fragmentary 
evidence  available.  The  principal  observations  have  been  briefly  summarized  and 
it  is  clear  that  they  are  contradictory  in  some  instances.  Before  an  attempt  is  made 
to  summarize  the  experimental  evidence.  Section  2  of  this  paper  will  describe  some 
new  observations  relative  to  the  pressure  distribution.  Section  4  is  an  att*'-  lpt  to 
assess  the  effect  of  the  observed  pressure  gradients  onvarious  types  of  measurements. 
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2.  EXPERIMENTAL  OBSERVATIONS 


A  dynamic  type  of  opposed  anvil  apparatus  has  been  in  use  at  AFCRL  for 
several  years.  The  design  of  the  apparatus  incorporates  several  novel  features  which 
have  been  described  in  the  literature  by  Riecker  (1964).  The  shear  press  is  used 
in  a  comprehensive  program  of  shear  strength  measurements  on  a  variety  of  rock 
and  mineral  materials  under  geophysically  realistic  conditions  of  temperature  and 
pressure  (Riecker  and  Seifert,  1964a,  1964b), 

A  series  of  static  compression  tests  were  conducted  on  compressed  mineral 
wafers  approximately  0.020  in.  thick  and  1/4  in.  or  1/2  in.  ?n  diameter.  The 
wafer^  vere  compacted  from  the  mineral  olivine  after  it  had  been  ground  to  pass 
a  400- mesh  screen.  The  maximum  and  median  grain  sizes  were  0.037  mm  and 
0.0043  mm,  respectively. 

Bismuth  and  thallium  wires  0.  010  in.  in  diameter  and  0  014  in.  long  were 
placed  vertically  in  the  prepared  wafers  at  various  dis\  rices  from  the  sample  cen¬ 
ter.  Thin  platinum  sheets  were  sandwiched  above  and  below  the  mineral  wafer 
to  complete  the  electrical  calibration  cell.  Resist,  nee  discontinuities  indicate 
the  bismuth  I- II  and  thallium  II- III  phase  transformations  at  25.  4  and  36.  7  '’-ilobars, 
respectively.  The  nominal  applied  pressure  required  to  produce  the  known  tran¬ 
sitions  in  the  calibrant  wires  was  recorded  as  a  function  of  the  radial  position  of 
the  wire.  Observed  pressures  were  divided  by  the  appropriate  transition  pressure 
in  order  to  form  a  dimensionless  parameter.  This  parameter  and  its  reciprocal 
are  plotted  in  Figure  1  against  the  sample  position  with  the  latter  expressed  as  a 
fraction  oil  the  anvil  radius.  It  is  evident  that  pressure  gradients  were  present. 

It  should  be  noted  that  when  expressed  in  this  normalized  fashion  the  data  ob¬ 
tained  on  the  bismuth  and  thallium  fall  on  a  single  line.  This  implies  that  the 
pressure  at  a  given  radial  position  is  proportional  to  the  nominal  applied  pressure. 

A  similar  observation  was  made  by  Myers  et  al  over  a  wider  pressure  range.  The 
connection  between  Figure  1  and  the  actual  pressure  gradients  will  be  developed 
in  the  next  section,  but  first  some  qualitative  observations  need  to  be  added. 

In  the  static  compression  experiments  the  presence  of  large  pressure  gradients 
was  also  revealed  by  button-like  central  indentations  produced  in  the  anvils,  similar 
to  the  effects  described  by  Bridgman(1937).  This  manifestation  of  the  pressure 
gradient  has  only  been  observed  in  static  tests;  in  shear  measurements  this  sort  of 
permanent  deformation  of  the  anvils  develops  very  rarely.  No  indentation  occurs 
in  static  compression  tests  when  the  sample  wafer  thickness  is  less  than  0.  018  in. 
for  l/2-in.-diam  anvils  or  0.014  in.  for  1/4-in. -diam  anvils.  This  observation 
agrees  qualitatively  with  that  made  by  Myers  et  a3;  me  pressure  multiplication 
effect  is  greatest  for  sample  geometries  with  small  diameter-to-thickness  ratios. 
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Figure  1  *  A  Normalized  Plot  of  the  Experimental  Pressures  at  which  Known 
Phase  Transitions  Were  Observed  to  Occur  as  a  Function  of  the  Radial  Posi¬ 
tion  of  the  Calibrant  Wire 


Qualitative  observations  made  with  the  AFCRL  apparatus  suggest  that  under 
dynamic  conditions  the  superimposed  shear  strain  greatly  reduces  the  pressure 
gradient  which  might  otherwise  be  present  in  the  sample.  In  addition,  high  temper¬ 
ature  usually  sign  ficantly  reduces  the  gradient  magnitude.  For  example,  in  static 
compression  tests  on  enstatite,  anvil  deformation  was  observed  at  room  temper¬ 
ature,  but  at  temperatures  of  300 °C  and  above,  significant  deformation  did  not  occur, 
even  at  pressures  up  to  75  kilobars.  WTien  test  wafers  are  first  heated  at  zero 
pressure,  sheared  also  at  zero  pressure,  and  then  raised  to  the  test  pressure,  the 
pressure  gradients  appear  to  be  minimal.  This  sequence  permits  the  efficient 
extrusion  of  excess  sample  material  early  in  the  pressurization  and  allows  a  wafer 
to  develop  an  equilibrium  thickness  of  from  six  to  eight  nils.  With  this  technique, 
any  rea  'onable  diameter- to-thickness  ratio  may  be  used  with  satisfactory  results. 
Samples  which  initially  ranged  in  thickness  from  0.010  in.  to  0.040  in.  yielded 
equivalent  shear  strength  or  friction  measurements  within  ±2  percent  when 
sheared  on  either  1/4- in.  or  l/2-in.-diam  anvils.  Unfortunately  similar  pro¬ 
cedures  cannot  be  used  to  eliminate  the  pressure  gradients  in  static  compres¬ 
sion  tests.  In  the  static  tests  the  actual  sample  configuration  must  be  evaluated 
experimentally  to  ensure  that  pressure  gradients  are  minimized. 
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On  the  basis  of  the  limited  evidence  available,  the  following  generalizations 
are  made: 

1.  Using  sample  materials  and  geometries  commonly  employed,  large  pres¬ 
sure  gradients  may  occur  in  static  anvil  devices. 

2.  The  maximum  pressure  usually  occurs  at  the  center  of  the  anvils  and  may 
be  as  much  as  2.  5  times  the  nominal  applied  pressure. 

3.  The  magnitude  of  the  pressure  multiplication  effect  is  markedly  dependent 
on  the  sample  geometry,  and  even  a  pressure  deficiency  at  the  anvil  center  may 
occur  in  thin  sample  wafers.  In  the  latter  case  the  maximum  pressure  will  occur 
in  an  annular  ring  concentric  with  the  axis  of  symmetry  as  predicted  by  Bridgman. 

4.  Frequently  the  pressure  at  a  given  radial  position  is  proportional  to  the 
nominal  applied  pressure  over  a  wide  range  of  loads. 

5.  A  transition  region  occurs  around  the  periphery  of  a  compressed  sample 
wafer  in  which  very  large  pressure  gradients  are  to  be  expected  because  this 
region  does  not  support  its  share  of  the  applied  load.  The  details  of  the  pressure 
distribution  in  this  region  are  co  upletely  unknown  at  present. 

6.  Direct  evidence  concerning  the  pressure  gradients  in  dynamic  anvils  is 
practically  nonexistent  because  of  severe  experimental  difficulties.  However  in¬ 
direct  eviuence  cited  previously  suggests  that  pressure  gradients  occurring  under 
dynamic  conditions  are  much  smaller  than  in  static  experiments. 

EMPIRICAL  ANALYSIS 

Several  theoretical  attempts  have  been  made  to  determine  the  pressure  gra¬ 
dients  existing  in  opposed  anvils  under  static  conditions  (Jamieson  and  Lawson, 

1962;  Jackson  and  Waxman,  1962;  Jackson  and  Davis,  to  be  published).  These  have 
not  yet  proved  to  be  of  great  value.  The  present  writers  approach  the  problem 
empirically,  realizing  that  although  the  results  will  be  less  satisfying  esthetically, 
the  empirical  approach  may  be  better  suited  to  their  purpose. 

The  analysis  proceeds  in  two  steps.  The  first  is  to  find  an  empirical  expression 
representing  the  pressure  distribution  between  compressed  anvils  which  is  simple 
enough  to  be  analytically  useful,  and  complex  enough  to  be  compatible  with  the  ex¬ 
perimental  evidence.  The  second  step  is  to  perform  calculations  using  the  em¬ 
pirical  equation  for  the  pressure  distribution  in  order  to  determine  the  magnitude 
of  the  errors  introduced  into  various  types  of  anvil  experiments  by  the  pressure 
gradients.  Shear  strength  l  leasurements  are  the  primary  concern  and  the  aim  is 
to  determine  the  order  of  magnitude  of  the  errors  which  can  arise;  for  example, 
whether  the  computed  shear  strengths  can  be  considered  accurate  to  5  or  50  per¬ 
cent.  Thus  certain  assumptions  which  are  obviously  over-simplifications  will  be 
made  in  the  interest  of  analytical  simplicity. 
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The  simplest  reasonable  assumption  which  can  be  made  is  that  in  a  compressed 
sample  wafer  the  pressure  varies  linearly  with  radial  position.  Such  a  relationship 
was  found  valid  experimentally  by  Montgomery  et  al,  although  their  evidence  is 
not  necessarily  unambiguous.  Following  Montgomery  et  A,  it  will  also  be  assumed 
that  the  pressure  drops  in  a  linear  manner  at  the  edge  of  the  sample  wafer  to  some 
small  value  comparable  to  the  shear  strength  of  the  material.  The  pressure  dis¬ 
tribution  can  be  represented  analytically  as 


P(p) 

Pc<l-«p> 

«  ^  P  Pe 

(la) 

Pip) 

pc  \T^r7(1-'» 

Pe<P^  1. 

(lb) 

where  p  =  — ,  that  is,  the  radial  position  is  expressed  as  a  fraction  of  V  e  anvil 
a 

radius,  a;  Pc  is  the  pressure  at.  p  0;  and  a  is  a  constant  which  determines  the 
pressure  gradient.  In  writing  Eq.  (lb)  it  is  assumed  that  the  shear  strength  of  the 
material  is  negligible  compared  to  the  pressures  being  applied,  so  that  the  pressure 
at  p  1  is  taken  to  be  zero.  This  is  a  reasonable  approximation  for  use  in  the 
high-pressure  realm  of  interest.  The  assumed  pressure  distribution  is  shown 
schematically  in  Figure  2. 

If  the  nominal  applied  pressure  is  determined  experimentally,  as  is  generally 
the  case,  the  central  pressure.  Pc,  can  be  expressed  as  a  function  of  a  and  p^  by 
integrating  the  pressure  distribution  over  the  anvil  area  and  equating  this  quantity 
to  the  total  applied  force.  Thus: 

i  ! 

Pa  T  f  P(P)  2npdp  .  (2) 

7T 

7r  a  0 

Using  Eq.  ( 1  )  it  is  readily  found  that 

Pc  -  f<  a,  flc>  P  (3a) 


where 

U°’  Pe>  l  f(i-a)pe(l+pe)  ‘ 

Then  fc-q.  (1  )  can  be  rewritten  as 

P(p)  =  fv«,  pj  P{1-  ap)  0  <p<pe 


(3b) 


(4a) 
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Figure  2 .  The  Sample  Pressure  as  a  Function  of  Radial  Position  as 
Calculated  from  an  Empirical  Equation  Derived  in  the  Text 


P(p)  -  g(o,  Pe)  P(1  -p)  Pe<P<  1 


(4b) 


where 


/  1  -  op  \ 

g(«.Pe>=  f("' pe)  •  (5) 

\  ^e  / 

Although  the  experimental  evidenc'  is  scant,  an  attempt  will  be  made  to  determine 
the  validity  of  the  simple  model  being  considered  by  comparing  Eq.  (4)  with  actual 
observations. 

In  experiments  where  calihrant  wires  are  located  at  various  radial  positions 
Eq,  (4)  predicts  a  linear  relationship  between  the  ^atio  of  transition  pressure  to 
nominal  applied  pressure  and  the  radial  position  over  the  central  portion  of  the 
anvils,  that  is, 

pt 

~  =  f<«.  p  )  (1  -atp)  p^  po .  (6) 

p  “  e 
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From  Figure  1  it  appears  that  the  data  obtained  with  the  AFCHL  press  approxi¬ 
mately  meet  this  requirement.  However,  there  is  some  indication  that  the  re¬ 
ciprocal  function  P/P^  results  in  a  better  straight-line  fit,  which  would  have  bv.en 
predicted  if  the  pressure  distribution  had  been  assumed  to  be  given  by 


p(p>  -  pc 


T+Tp 


pe 


(7a) 


(  1  -  p) 

P(P)  =  P c  (  1  +  0  pe  )  (  1  -  pe  )  Pe<P^l. 


(7b) 


Using  £q.  (7  )  and  balancing  forces,  IJc  can  be  eliminated  as  above.  Thus 


P(p) 


P(p) 


h(/3,  pe  )  (  1  ”  P  pl 


P( 1 -p) 


p<  Pe 

pe<  p<  1 


(8a) 


(8b) 


where 


m  &  Pe  i  ^  [0  pe 


In  ( 1  +  0  pe )] 


(9) 


and 


i  ( P.  Pe>  s  h  (0,  pe  )  (  1  +  j3  pe)  (  1  -pe )  . 

Then  the  equation  analogous  to  Eq.  (6)  is: 

h(fi,  p^)  (  1  +/3  p)  p<  pe  . 


(10) 


(ID 


Therefore  over  the  central  porfior.  of  the  anvils  the  pressure  distribution  assumed 
inEq.  (7)  results  in  a  linear  relationship  between  the  ratio  F/P  and  p  as  expressed 
in  Eq.  ( 1 1). 
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The  data  appear  to  fit  Eq.  (11)  slightly  better  than  Eq.  (S),  although  the  scatter 
and  reproducibility  of  the  data  are  such  that  a  significant  choice  cannot  be  made 
between  these  alternative  formulations.  The  measurements  maJe  by  Myers  et  al 
were  made  at  only  a  few  values  of  p  and  show  more  scatter  than  the  present  results 
so  that  their  data  also  do  not  permit  a  decisive  choice.  The  measurements  by 
Montgomery  et  al,  however,  are  much  more  accurate.  In  their  initial  measure¬ 
ments  using  bismuth  calibrant  wires  it  was  shown  that  the  aata  fit  an  equation  of  the 
form  of  Eq.  (11)  very  well.  In  the  latter  measurements  using  manganic  wires  it 
was  shown  that  Eq.  (6)  was  applicable.  This  contradiction  is  more  apparent  than 
real,  however,  because  when  0  is  small,  Eq.(ll)  is  equivalent  to  Eq.  (6).  Thus  it  is 
not  possible  with  the  experimental  data  available  to  make  a  clear  choice  between 
the  alternative  pressure  distributions  suggested.  Clearly,  more  detailed  data  are 
needed.  However  there  is  no  need  at  present  to  assume  a  more  complex  pressure 
distribution  than  those  that  have  been  considered.  Accordingly,  the  linear  pressure 
gradient  given  byEq.  (4)  will  be  accepted  tentatively  as  valid  for  the  purpose  of  con¬ 
tinuing  calculations. 

It  should  be  noted  that  both  of  the  pressure  distributions  suggested  are  slightly 

unrealistic  in  that  they  predict  a  cusp  at  the  center  of  the  anvils.  Physically  it  is 
dP 

expected  that  ^  -  0  at  this  point.  Most  of  the  theoretical  developments  have  also 
suffered  from  this  difficulty,  which,  however,  is  not  serious  because  only  a  very 
smail  fraction  of  the  sample  area  at  the  center  of  the  anvils  is  uncertain.  It  might 
be  added  that  the  pressure  distribution  given  byEq.  (8)  qualitatively  resembles  the 
theoretical  predictions  (Kiecker  and  Seifert,  1964a). 

The  values  of  the  constants  a  and  pe  may  be  calculated  from  the  slope  and 
intercept  of  the  straight  line  fitted  to  the  data  in  Figure  1  according  to  Eq.  (6). 

These  values  were  found  to  be  about  0.  46  and  0.  84  respectively,  on  the  basis  of  a 
least- squares  fit  from  wh*ch  the  uncertainties  were  estimated  to  be  10  percent 
and  5  percent,  respectively.  The  radial  pressure  gradient  existing  in  the  sample 
wafer  can  be  determii  ed  from  Eq.  (4).  Thus 


dP  _  fa_  p 
dr  a 


0<  p<  pe 


and 


(12a) 


dP 

dr 


-El  p 

a 


pg<  p<  i . 


(12b) 
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Over  the  central  portion  of  the  anvils,  the  pressure  gradient  in  the  AFCRL 
press  is  found  to  be  -3.  4P  in  the  pressure  range  investigated.  Thus  at  a  nominal 
applied  pressure  of  30  kilobars  a  gradient  of  about  100  kilobars  per  inch  would  be 
present.  The  gradient  in  the  annular  ring  aro  nd  the  periphery  of  the  sample  is 
much  larger,  amounting  to  about  -39P.  These  results  may  be  compared  with  those 
of  Montgomery  et  al,  where  the  central  gradient  was  found  to  be  about  +0.  47P  and 
the  edb?  gradient  amounted  to  about  -33P.  The  small  central  gradient  they  ob¬ 
served  no  doubt  results  from  the  Ag  Cl  composition  of  the  central  disc  of  the  sample 
cell,  which  is  a  much  more  ductile  material  than  the  olivine  used  in  the  AFCRL 
tests  Similarly,  the  relatively  close  agreement  on  the  edge  gradient  is  probably 
connected  with  the  fact  that  the  pyrophyllite  retaining  rings  used  by  Montgomery 
et  al  have  mechanical  properties  more  co  iparable  to  olivine.  Montgomery  et  al 
also  cited  evidence  which  indicated  that  the  pressure  started  dropping  sharply 
about  0.030  in.  in  from  the  edge  of  the  sample  cell,  that  is,  the*r  results  indicate 
pe  =  0.  88  which  is  in  reasonable  agreement  with  the  AFCRL  measurements.  In 
addition  the  compression  tests  on  pyrophyllite  wafers  by  Montgomery  and  Vaisnys 
indicated  that  an  annular  ring  around  the  periphery  of  the  sample  supported  only  a 
small  fraction  of  the  applied  load.  For  a  wafe  of  the  same  initial  thickness  and 
diameter  as  those  used  in  the  AFCRL  experiments  their  data  indicate  that  29.  6  per¬ 
cent  of  the  sample  area  was  non- load- supporting,  that  is,  Pe  was  about  0.  84,  in 
satisfactory  agreement  with  the  other  observations  made  above. 

These  numerical  comparisons  indicate  that  the  pressure  distribution  given 
byEq,  (4)  is  n^t  seriously  in  error.  This  pressure  distribution  is  shown  in  Figure 
2,  which  uses  the  values  of  a  and  pe  deduced  from  the  AFCRL  data,  assuming  a 
nominal  applied  load  of  25.4  or  36.7  kilcbars.  Note  that  a  pressure  multiplication 
effect  of  about  1.  6  was  observed. 

The  ioregoing  considerations  lead  to  the  conclusion  that  the  existence  of  a 
large  pressure  multiplication  effect,  as  indicated  by  data  such  as  that  shown  in 
Figure  1,  does  not  necessarily  imply  the  existence  of  a  large  pressure  gradient 
over  the  central  portion  of  the  sample  cell.  It  may  only  mean  that  there  is  an 
annular  region  a  j  the  periphery  of  the  sample  which  is  not  supporting  its  share 
of  the  applied  load.  This  analysis  and  the  available  data  do  indicate  that  the  magni¬ 
tude  of  the  pressure  gradient  over  the  central  portion  of  the  sample  increases  in 
proportion  to  the  total  applied  load,  and  that  the  actual  gradient  present  in  a  given 
situation  will  strongly  depend  on  the  composition  and  dimensions  of  the  sample  cell 
as  well  as  on  the  total  applied  load.  It  is  also  clear  that  very  large  gradients  are 
to  be  expected  at  the  sample  edge  where  a  ring  comprising  as  much  as  30  percent 
of  the  totai  anvil  area  may  be  supporting  only  a  small  fraction  of  the  applied  load. 
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3.1  Static  \n\il  Experiments 


The  influence  of  a  pressure  distribution  ruch  as  that  given  byEq.  (4)  on  various 
types  f  experiments  will  now  be  considered. 

3.  1.  1  PHASE  TRANSITIONS  II\  HOMOGENEOUS  DISCS 

First  consider  a  case  where  the  central  portion  of  the  sample  wafer  is  a  homo¬ 
geneous  material  which  undergoes  a  phase  transition  upon  the  application  of  high 
pressure.  If  a  is  positive,  as  is  commonly  the  case,  the  pressure  is  greatest  at 
the  center  of  the  sample  wafer  and  the  phase  transition  begins  at  the  center  when 
the  nominal  applied  pressure  reaches  a  value  given  by 


where  is  the  true  pressure  at  which  the  phase  transition  occurs. 

As  the  pressure  is  increased  the  transformed  region  continues  to  increase  in 
diameter  until  it  extends  over  the  entire  area  bounded  by p^.  The  nominal  applied 
load  then  is  given  by 


*  f(l-ape)  * 


(14) 


Further  increasing  the  pressure  will  result  in  only  a  very  slow  increase  in  the  area 
of  the  transformed  material  because  of  the  very  sharp  pressure  gradients  existing 
near  the  edge  of  the  sample  cell.  Thus  the  width  cf  the  phase  transition,  as  re¬ 
flected  in  the  range  of  nominal  applied  pressures,  required  to  drive  it  to  completion 
can  be  determined  from  the  difference  of  Eqs.(13)  and  (14).  Thus 


AP  = 


(15) 


Note  that  this  transition  interval  is  due  only  to  the  pressure  gradients;  the  actual 
transition  is  assumed  to  be  perfectly  sharp.  In  addition,  the  transition  pressure 
determined  by  this  hypothetical  experiment  can  be  determined  by  averaging  Eqs. 
(13)  and  (14).  Thus 


2 -  «Pe 


2f 


(if.) 


JCV 
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where  is  the  transition  pressure  as  determined  in  t\e  hypothetical  experiment. 

Using  the  values  of  Qt,  p^,  and  f  determined  in  the  AFCRL  experiments  as 
illustrative  of  a  typical  case,  it  is  found  that 


p.  - 

1 

0. 61  P  , 

(17a) 

V 

0.  99  P  , 

(17b) 

AP  = 

0. 38  Pt , 

(17c) 

V 

0.  80  P 

( 17d) 

Thus  the  transition  begins  when  the  nominal  applied  pressure  is  only  61  percent  of 
the  true  transition  pressure  and  proceeds  until  it  reaches  99  percent  of  the  true 
pressure.  The  transition  would  appear  sluggish,  since  the  wHth  would  be  about 
38  percent  of  the  transition  pressure,  and  the  average  transition  pressure  deduced 
from  this  hypothetical  experiment  only,  would  be  80  percent  of  the  true  value. 

Although  this  numerical  example  cannot  be  taken  too  seriously,  it  illustrates 
clearly  the  nature  of  the  problems  which  can  be  encountered  in  the  investigation  of 
opposed  anvils.  The  numerical  values  used  to  characterize  the  pressure  distribu¬ 
tion  were  shown  previously  to  be  generally  compatible  with  tnose  reported  by  other 
investigators.  It  is  of  interest  to  determine  how  these  results  would  be  affected 
by  reasonable  variations  in  the  parameters  0 1  and  p^. 

It  is  clear  that  if  a  is  zero,  the  width  of  the  transition  will  also  be  zero,  but 
that  the  observed  transition  pressure  will  still  be  less  than  the  true  value  because 
f  has  a  value  greater  than  1.0  for  any  p^  less  than  1.  0.  The  detrimental  effect  of 
the  nor- load-bearing  circumferential  region  is  very  evident..  All  of  the  experi¬ 
mental  evidence  available  suggests  that  very  little  canbe  done  experimentally  to  in¬ 
crease  p^  above  0.9,  although  this  factor  depends  somewhat  on  the  sample  material  and 
dimensions.  It  is  therefore  clear  that  the  nominal  applied  load,  which  is  the  quantity 
usually  observed  experimentally,  mayb°  considered  a  parameter  which  mustbe  cali¬ 
brated  using  known  phase  transitions  as  fixed  points,  or  by  some  other  means.  When  a 
is  negative,  as  it  was  in  the  experiments  of  Montgomery  et  al,  the  error  introducedby 
the  negative  pressure  gradient  over  the  central  portion  of  the  sample  may  compen¬ 
sate  part' ally  for  the  error  introduced  by  the  non- load-bearing  peripheral  region. 

In  a  more  common  experimental  arrangement  the  sample  undergoing  the  phase 
transition  is  surrounded  by  a  retaining  ring  of  some  other  material.  In  this  case, 
Pe  should  be  replaced  b^  pf,  that  is,  by  a  value  indicative  of  the  boundary  between 
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the  sample  and  the  retaining  ring.  Inasmuch  as  is  usually  less  than  th  width 
of  the  transition  would  appear  smaller,  but  the  average  transition  pressure  would 
be  even  more  in  error.  In  this  modified  example,  the  value  of  f  is  not  changed  as 
long  as  it  is  assumed  that  the  pressure  gradient  remains  the  same. 

Another  factor  which  might  affect  the  preceding  analysis  is  the  influence  of  any 
volume  change  associated  with  the  phase  transition.  As  the  phase  transition  pro¬ 
ceeds  under  increasing  pressures  the  transformed  material  generally  undergoes  a 
reduction  in  volume,  and  hence  it  tends  to  withdraw  support  from  the  anvil.  This 
causes  an  increase  in  the  pressure  applied  to  the  low-pressure  polymorph,  and  tends 
to  drive  the  transformation  to  comnletion  somewhat  faster  th  in  was  predicted  by 
the  simple  model  considered.  Although  the  model  could  be  made  more  sophisti¬ 
cated  by  considering  volume  effects,  the  numerical  data  available  and  the  accuracy 
of  the  calculations  that  can  be  made  do  not  justify  the  additional  effort  at  present. 

3  1.2  ELECTRICAL  RESISTANCE  DISCONTINUITIES 

As  a  second  example,  consider  the  case  in  which  a  resistance  element  in  the 
form  of  a  wire  or  narrow  ribbon  is  placed  across  the  diameter  of  the  anvil  face. 

It  is  assumed  that  the  wire  is  embedded  in  silver  chloride  or  some  other  suitable 
pressure-transmitting  medium,  and  that  the  pressure  dependence  of  the  resistance 
of  the  sample  is  observed.  Consider  in  particular  that  the  materia  undergoes  a 
phase  transition  with  which  there  is  associated  a  discontinuity  i;?  electrical  resis¬ 
tivity.  This  is  the  experimental  arrangement  conmonl)  used  by  Bridgman  and 
others. 

The  analysis  given  previously  is  directly  applicable  to  this  situation.  It  is 
highly  unlikely  that  phase  transitions  could  be  located  with  precision  from  measure¬ 
ments  of  nominal  applied  pressure  at  which  resistance  discontinuities  are  observed. 
Thus  it  is  not  surprising  that  Bridgman's  anvil  measurements  disagreed  seriously 
with  his  piston  and  cylinder  measurements  where  more  nearly  hydrostatic  conditions 
prevailed.  The  pressure  distribution  in  an  anvil  geometry  approximating  that  used 
by  Bridgman  was  measured  by  Montgomery  et  al  and  was  discussed  earlier. 
Montgomery  et  al  found  that  th«*  discrepancies  in  Bridgman’s  measurements  could 
easily  be  accounted  for  by  the  observed  pressure  gradients. 

3.2  Dynamic  Anvil  Experiments 


The  preceding  discussion  referred  speeificallj  to  static  anvils.  When  dynamic 
anvils  are  considered  the  problems  become  more  complex.  For  example,  the 
occurrence  of  a  phase  transition  might  be  as  described  previoi^sly  in  some  simple 
cases;  however,  shear  strain  frequently  facilitates  phase  transitions  and  adds  a 
complicating  factor.  In  fact  the  effect  of  shear  strain  on  reaction  rate  provides 
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the  primary  motivation  (or  performing  phase  diagram  studies  and  synthesis  work 
on  dynamic  anvils. 

An  important  application  of  dyna.o'c  anvils  is  the  measurement  of  the  shear 
strength  of  materials  at  high  pressures.  The  AFCRL  experiments  have  been  dis¬ 
cussed  in  detail  elsewhere  (Riecker  and  Seifert,  1964a,  1964b);  however  the  essen¬ 
tial  features  need  to  be  described  here  briefly. 

The  sample  consists  of  a  homogeneous  disc  compressed  between  the  anvils. 

One  anvil  is  rotated  at  a  constant  rate  by  a  suitable  drive  mechanism  while  the  other 
is  held  stationary.  The  torque  required  to  shear  the  sample  is  measured  as  a 
function  of  the  applied  load,  temperature,  and  rate  of  strain.  From  these  data  and 
the  known  geometrical  factors,  the  shear  strength  of  the  sample  material  may  be 
determined  as  a  function  of  pressure,  temperature,  and  rate  of  strain. 

Experimentally  it  is  found  that  at  low  pressrres  the  sample  does  not  shear  but 
rather  that  surface  slippage  occurs  at  the  interface  between  the  sample  and  anvil. 

As  the  pressure  is  increased,  a  point  is  reached  at  which  th»  surficial  friction 
stresses  exceed  the  shear  strength  of  the  sample,  and  internal  shearing  begins.  The 
pressure  required  to  produce  the  transition  from  slippage  to  shearing  varies  greatly 
with  the  properties  of  the  sample  material;  for  the  hard  rock  and  mineral  materials 
with  which  the  authors  are  particularly  concerned,  shearing  does  not  begin  until 
pressures  of  about  30  kilobars  are  reached.  A  typical  plot  of  turning  moir.er  t,  M, 
versus  nominal  applied  pi  ssure,  P*  is  shown  in  Figure  3. 

3.  2.  1  SURFICIAL  FRICTION  MEASUREMENTS 

Consider  first  the  low-pressure  region  in  which  surficial  slippage  occurs. 

Using  the  pressure  distribution  given  by  eq  (4)  and  assuming  the  coefficient  of 
friction  between  the  sample  material  and  a to  be  independent  of  pressure,  the 
turning  moment,  M,  can  be  found  from 


M  2tt  a'^  jj.  f  P(p)p^dp 
0 


(18a) 


or 


(18b) 
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Figure  3.  A  Typical  Plot  of  Turning  Moment  versus  Nominal  Applied  Pressure 
as  Observed  in  Shearing  Experiments.  The  knee  in  the  curve  separates  the  low 
pressure  realm  in  which  surficial  slippage  occurs  from  the  high  pressure  realm 
in  which  internal  shearing  of  the  sample  occurs 


This  reduces  to 

MP  jfPe3  d-|«Pe)  +g[l-Ps3<1-fPe)]j  '  (19) 

where  the  symbols  have  the  same  meanings  used  previously.  Thus  the  friction 
coefficient  can  be  calculated  from 


U  = 


3Q(or.p  ) 


(20) 


where  Q(  Pe )  defined  to  be  the  reciprocal  of  the  expression  inbxackets  in 
Eq.  (19). 

Equation  (20)  is  well  established  experimentally.  Measurer^  nts  made  on  min¬ 
eral  materials  show  ihe  predicted  linear  increase  in  turning  moment  with  pressure. 
In  addition  the  numerical  values  for  the  friction  coefficient  agree  with  published 
data.  Bridgman  reported  similar  observations  on  a  wide  range  of  sample  materials. 


18 


The  quantity  Q(ct,pe)  is  a  correction  factor  which  would  be  l.U  in  the  idea) 
case  of  a  completely  uniform  pressure,  that  is,  when  a  =  0  and  =  1.  0.  Tn  any 
real  experiment  Q  will  deviate  somewhat  from  unity,  the  extent  of  this  deviation  is 
the  subject  of  interest.  Numerical  values  of  Q  are  shown  graphically  in  Figure  4 
for  a  range  of  values  of  a  and  Using  the  values  of  these  parameters  derived 
from  the  static  AFCRL  measurements,  it  is  found  that  Q  ~  1.  10,  which  means  that 
friction  coefficients  calculated  assuming  Q  =  1.  0  will  be  too  small  by  cbout  10  per¬ 
cent.  However,  this  is  a  pessimistic  estimate.  The  experimental  observations 
discussed  in  Section  2  indicate  that  a  is  smaller  and  p^  larger  under  dynamic  con¬ 
ditions.  This  probability  is  further  supported  by  analysis  of  shear  measurements 
given  later.  Thus  it  is  reasonable  to  conclude  that  friction  measurements  made 
in  this  manner  will  not  be  in  error  by  more  than  a  few  percent. 

3.  2.  2  SHEAR  STRENGTH  MEASUREMENTS 

In  the  high-pressure  region  surficial  slippage  ceases  and  the  sample  wafer 
e  periences  internal  shearing.  The  shea  strength  of  the  material  may  be  expected 
to  be  dependent  on  the  pressure,  temperature,  and  rate  of  strain  in  the  general  case. 
In  order  to  analyze  the  problem  several  simpl  'ying  assumptions  must  be  made. 

Thus  it  is  assumed  (1)  that  the  shear  strength  of  the  material  increases  linearly  wuh 
the  applied  pressure  (2)  that  the  temperature  is  uniform  throughout  the  sample 
wafer,  and  (3)  that  the  shear  strength  is  independent  of  the  rate  of  strain  of  the 
sample  material. 

The  first  assumption  is  generally  compatible  with  the  results  obtained  here 
and  with  numerous  observations  reported  by  Bridgman.  It  cannot  be  expected  to 
hold  in  the  vicinity  of  phase  transitions,  however.  The  second  assumption  is 
reasonable  in  view  of  the  small  thickness  of  the  samples  and  the  relatively  high 
thermal  conductivity  and  large  mass  of  typical  anvils.  The  third  aosumption  is  not 
generally  valid.  It  was  shown  by  Bridgman  (193?)  that  the  shear  strength  of  some 
materials  was  independent  of  rate  of  strain  if  the  experimental  temperature  was 
much  lower  th  the  melting  point  of  the  material.  This  condition  is  generally 
satisfied  in  the  AFCRL  experiments  and  only  small  rate  effects  have  been  observed 
over  several  orders  of  magnitude  of  anvil  speed.  The  third  assumption  must  be 
used  witu  caution,  however,  because  the  rate  of  strain  of  the  sample  wafer  varies 
from  essentially  zero  at  the  axis  of  rotation  to  some  large  value  at  the  periphery, 
depending  on  the  rate  of  anvil  rotation. 

The  transition  from  surficial  slippage  to  internal  shearing  does  not  occur 
abruptly.  Due  to  the  pressure  gradients  present  in  the  sample,  a  transition  interval 
is  encountered  which  results  in  the  knee  shown  in  Figure  3;  the  larger  the  gradients 
the  wider  the  transition  region. 
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Figure  4.  The  Correction  Factor  Q(ottpQ)  as  a.  Function  of  the  Parameters 
ot  andpe<  This  factor  is  defined  in  the  text,  Eq.  (20) 


If  it  is  assumed  that  the  pressure  distribution  given  by  Eq.  (4)  is  applicable, 
then  the  internal  shearing  will  begin  at  the  center  of  the  sample  when  pP^-T 
where  r  is  the  shear  strength  of  the  material.  This  is  assumed  to  be  linearly  re¬ 
lated  vO  pressure  as 

T=S^P,  (21) 

where  S  is  the  shear  strength  oi  the  material  in  the  absence  of  hydrostatic  pressure 
and  Q  is  a  constant.  S  and  0  may  be  dependent  on  temperature  but  are  assumed  to 
be  independent  of  rate  of  strain. 

Under  these  conditions  it  is  found  tnat  internal  shearing  will  start  at  the  center 
when  the  nominal  applied  pressure  reaches  the  value 

=  (  jrjTf  •  (22> 

As  the  pressure  is  increased  the  diameter  of  the  sheared  portion  will  increase 
rapidly  until  pP  (pg)  =  r,  that  is,  until  the  radius  of  the  sheared  region  reaches  /a ) 
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Further  increases  in  pressuie  will  result  in  a  more  gradual  increase  in  the 
shearing  area  because  of  the  sharp  drop  in  pressure  near  the  sample  edge.  Thus 
the  transition  to  shearing  is  nearly  com  'ete  when 


Pc  = 


f  ( p  ii  )  (  1  -  ap  )  f 


(23) 


The  width  of  the  transition  is  readily  found  to  be 


or  o 


AP 


\li-0)  (1  -op  >f 


(24) 


and  the  average  pressure  at  which  the  transition  occurs  is  given  by 


2- 


P.  = 


N  (M-/m  \  1-flfft, 


(25) 


Thus  the  ratio  of  these  quantities  is 


AP 


2  ap 


"R  <«.p€>=  — 


N 


aP~ 


(26) 


The  ratio  R(Of,  )  is  shown  in  Figure  5  for  an  appropriate  range  of  the  param¬ 
eters.  Note  that  this  ratio  is  dependent  on  the  form  of  the  pressure  gradient,  but 
is  not  explicitly  dependent  on  the  properties  of  the  material  being  sheared.  Equation 
(26)  is  also  applicable  to  phase  transformations  provided  the  transformation  is  not 
influenced  shear  or  is  observed  under  static  conditions. 

Numerical  values  of  the  ratio  R(cr,Pe)  can  be  obtained  from  experimental  data 
which  will  permit  some  limits  to  be  placed  on  the  pressure  distribution  under  dv- 
nani  c  conditions.  For  example,  the  AFCRL  measurements  on  minerals  indicate 
that  AP  is  definitely  less  than  5  kilobars  while  P^T  is  about  30  kilobars.  Thus  R 
is  less  than  0.  17  in  these  experiments.  Reference  to  Figure  5  reveals  that  for 
p '■  in  the  range  0.9  to  1.0,  at  is  less  than  0. 17,  that  is,  much  less  than  observed  in  the 
static  experiments  onolivine.  Thus  the  width  and  location  of  the  knee  in  the  plot  of 
turning  moment  versus  nominal  applied  pressure  provides  an  internal  check  on  the 
pressure  gradients  present  in  the  shear  experiments  and  corrobeivtes  the  quali¬ 
tative  observations  made  earlier. 
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Figure  5.  The  ratio  H(o/t  p  )  as  a  Function  of  the  Parameters  a  and  p  .  This 

e  0 

factor  is  defined  in  the  text,  Eq.  (26) 


The  turning  moment  required  to  shear  the  sample  wafer  may  be  computed  from 
the  integral 


M  =  2k  a3  f  r(P)p2d p.  (27) 

0 

Considering  the  high-pressure  realm  where  the  radius  of  the  shearing  region  is 
Pe  or  greater,  and  using  the  usual  pressure  distribution,  this  integral  can  be  written 
as 


2k  a3 


[S  +  /3  fP  (  1  -  op)]  p2  dp 


_S  +  ft  gP(  1  -p)]  p2  dp 


+  /  pgP(  1  -  p)  p2  dp 


(28) 
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Here  p^  is  the  radius  marking  the  outer  boundary  of  the  shearing  region.  Although 
this  general  expression  can  be  readily  integrated,  it  is  more  convenient  for  numer¬ 
ical  calculations  to  consider  two  special  cases,  that  is.  when  Px  '  Pe  and  when 
p  -  1.0  The  former  case  corresponds  to  the  most  unfavorable  situation  which 
applies  above  the  "knee"  in  the  turning  moment  curve  (Figure  3).  In  this  case  only 
surfici;-T  slip  occurs  in  the  annular  region  between  P  rPe  and  p  1.0.  The  second 
case  to  be  considered  applies  in  the  very-high-pressure  realm  where  the  material 
in  the  annular  ring  is  also  being  sheared. 

In  the  first  case  the  turning  moment  is: 


M 


2tt  a 


Pe3  5  +  P  pe3  f  (  1  '  +  ^  [\  -  Pe 3  <  1  '  f  Pe )]  ( 


L 


while  in  the  latter  case  it  takes  the  simpler  form: 


(29) 


S  -t  3 


f  (  1  - 1  o  Pe  )  +  g 


(30) 


The  linear  relationship  between  M  and  P  predicted  by  Eqs.(29)  or  (30)  is  in 
agreement  with  measurements  made  on  mineral  materials  in  the  AFCRL  app  ratus 
as  well  as  with  many  measurements  ma^0  by  Bridgman  (1935,  1936,  1937).  The 
constants  S  ana  *3  can  be  determined  from  the  intercept  and  slope,  respectively,  of 
a  plot  of  turning  moment  versus  nominal  applied  pressure.  When  these  graphically 
determined  factors  are  inserted  into  Eq.  (21)  the  piessure  dependence  of  the  shear 
strength  of  the  material  under  test  is  found  from  the  following  two  expressions 
corresponding  to  the  two  special  cases  under  consideration: 


T  (P)  r, 

2?ra 


,M0  +  T  (<v,  pe,  ix  Z) 


dM 

dP 


or 


(31) 


T (P)  -2— 

2ia 


M, 


+  Q(  a,Pe) 


dM  P 
dP 


(32) 


In  these  expressions  is  the  intercept  found  by  extrapolation  as  indicated  in 

Figure  3,  QU'Y,  p  )  remains  as  previously  defined,  and  T(or,p  ,  pZ)  is  given  by 
e  e 


i 
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where 


(33) 


.  2?ra3  dlJ 
3  '  cfi\f 


(34) 


Consider  first  the  unfavorable  case  depicted  byEq,  (31).  Two  correction  factors, 

-  3 

p^  and  T,  are  encountered  which  are  normally  assumed  to  be  unity  in  the  calcu¬ 
lation  of  shear  strengths.  The  neglect  of  these  factors  can  lead  to  appreciable  errors, 

-  3 

depending  on  the  pressure  gradients  present.  For  example,  p^  tajtes  on  the  values 
1.37,  1.  17  and  1.00  when  p^  assumes  the  values  0.90,  0.95,  and  1.00,  .respectively. 
Thus  the  first  term  in  Eq.  (31)  may  be  in  error  by  as  much  as  37  percent  under 
unfavorable  conditions. 

The  function  T  pZ)  cannot  be  calculated  in  the  general  case  because  of 

the  presence  of  the  factor  Zp,  but  for  the  mineral  materials  of  particular  interest 
this  quantity  is  found  to  have  the  value  2.0.  For  this  special  case  the  function 
T ( a ,  p^ ,  pZ)  has  been  plotted  in  Figure  6  for  an  appropriate  range  of  the  parameters. 
Note  that  this  quantity  may  be  either  greater  or  less  than  unity.  For  a  reasonable 
range  of  parameters,  such  as  p^  >0.  9  anaa<0.  2,  it  is  clear  that  the  second  term 
inEq.  (31)  should  not  be  in  error  by  more  than  12  percent. 

These  considerations  lead  to  the  conclusion  that  shear  strengths  calculated  from 
Eq.(31)  with  the  assumption  that  p^  and  T  are  unity  will  be  too  small  by  as  much  as 
38  percent  at  very  low  pressures,  and  too  large  by  at  most  12  percent  at  very  high 
pressures.  At  intermediate  pressures  the  error  will  be  somewhere  between  these 
extremes  depending  on  the  relative  magnitude  of  the  two  terms  on  the  right  hand 
side  of  Eq.  (31).  It  is  also  apparent  that  the  error  arises  primarily  from  *he  fact 
that  p^  f  1.  Thus  an  annular  ring  of  material  around  the  periphery  of  a  sample 
which  is  not  sheared  can  be  ver^  detrimental. 

Consider  now  the  most  favorable  case  which  can  occur,  that  is,  when  px  ~  1.0 
and  the  entire  sample  is  sheared.  The  shear  strength  is  given  by  Eq.  (32)  and  there 
is  only  one  correction  factor,  Qfoqp^),  which  is  plotted  in  Figure  4.  For  reason¬ 
able  values  of  the  parameters,  that  is,pe>  0.  9  and  nr  <  0.  2,  ix  is  found  that  the  error 
in  the  second  term  of  Eq.  (32)  should  be  less  than  7  percent.  Thus  xhe  shear  strengths 
calculated  from  Eq.(32)  would  be  in  error  by  a  negligible  amount  at  very  low  pres¬ 
sures  while  they  would  be  too  small  by  as  much  as  7  percent  ax  very  high  pressures. 
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PARAMETER  a 


Figure  6.  The  Correction  Factor  T(  atp  ,  (jlZ)  at?  a  Function  of  the  Parameters 

e 

a  and  for  the  Special  Case  pZ  =  2.0.  These  factors  are  iefined  in  the  text, 
Eqs.  (33)  and  (34) 


On  the  basis  of  these  two  limiting  cases  it  is  reasonable  to  conclude  that  shear 
strength  measurements  made  in  dynamic  anvils  will  not  be  grossly  in  error,  even 
in  the  presence  of  rather  substantial  pressure  gradients.  In  the  measurements 
on  minerals  just  presented,  it  seems  probable  that  the  computed  shear  strengths 
are  not  in  error  by  more  than  10  percent.  However,  it  is  clear  that  the  peripheral 
region  around  the  edge  of  the  sample  disc  where  the  pressure  drops  i  apidly  has  a 
decisive  bearing  on  the  accuracy  of  these  calculations  as  well  as  on  those  discussed 
previously.  The  width  of  this  peripheral  region  as  characterized  by  p^  remains  a 
largely  unknown  but  crucial  factor,  especially  in  dynamic  anvils. 


4.  C.ONCU  SONS 

The  experimental  evidence  relative  to  the  pressure  gradients  in  opposed  anvil 
devices  was  reviewed  briefly  in  Section  1  and  some  new  observations  were  given  in 
Section  2.  The  data  were  found  to  be  generally  ^onsiste  \  A  compressed  sample 
may  be  divided  conveniently  into  two  major  r<.  ^  ns:  (1'  A  central  portion  com¬ 
prising  about  7  5  percent  of  the  anvil  area,  in  vhich  the  pressure  usually  decreases 


25 


toward  the  sample  edge,  although  a  negative  gradient  may  be  presert  under  some 
conditions,  for  example,  when  the  sample  cell  has  a  very  large  diameter- to- thi>k 
ness  ratio,  (2)  An  annular  ring  around  the  periphery  of  the  sample  comprising 
about  25  pe  cent  of  the  area,  where  the  pressure  drops  rapidly  to  some  small 
value  indicative  of  the  shear  strength  of  the  material. 

Section  3  gave  an  empirical  formula,  assuming  a  linear  gradient,  for  the  pres- 
sure  distribution  in  these  two  regions,  which  was  found  to  be  compatible  with  the 
experimental  observations.  This  pressure  distribution  was  used  to  estimate  the 
errors  which  might  be  introduced  into  varkr  ypes  of  anvil  experiments  by  the 
pressure  gradients.  It  was  found  that  appreciable  errors  oau  arise,  and  correction 
factors  were  given  graphically  for  a  variety  of  conditions.  Analysis  of  some  shear 
experiments  indicated  that  the  pressure  gradients  are  much  smaller  under  dynamic 
conditions  than  under  static  conditions. 

These  considerations,  although  generalized,  have  shown  that  in  anvil  experiments 
very  careful  consideration  must  be  given  to  the  effects  of  pressure  gradients.  In 
*  tatic  experiments  the  sample  geometry  may  be  adjusted  and  proper  calibration 
procedures  used  such  that  only  minor  errors  are  introduced.  Under  dynamic  con¬ 
ditions  very  little  can  be  done  to  modify  the  pressure  gradiems  and  it  is  fortunate 
that  the  shearing  of  the  sample  cell  tends  to  bring  about  a  more  uniform  pressure  distri¬ 
bution.  The  analysis  presented  has  shown  that  in  many  instances  the  peripheral  region, 
which  is  partially  non- load- supporting,  has  a  greater  detrimental  effect  on  anvil 
experiments  than  does  the  pressure  gradient  which  exists  over  the  central  portion 
of  the  sample  cell. 

The  empirical  analysis  in  Section  3  provides  sonic  insight  into  the  proolems 
which  can  arise  in  anvil  experiments  due  to  the  presence  of  pressure  gradients, 
and,  provided  ^ome  knowledge  of  the  pressure  gradients  is  available,  it  permits 
rough  computations  of  the  errors.  It  is  hoped  that  additional  systematic  investi¬ 
gations  ’ll  be  made  into  the  nature  and  magnitudes  of  pressure  gradients  so  that 
these  ideas  may  be  more  rigorously  applied  or  more  sophisticated  ones  developed. 
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